Bose-Einstein Condensate in Weak 3d Isotropic Speckle Disorder 
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The effect of a weak three-dimensional (3d) isotropic laser speckle disorder on various thermody- 
namic properties of a dilute Bose gas is considered at zero temperature. First, we summarize the 
derivation of the autocorrelation function of laser speckles in Id and 2d following the seminal work 
of Goodman. The goal of this discussion is to show that a Gaussian approximation of this func- 
tion, proposed in some recent papers, is inconsistent with the general background of laser speckle 
theory. In this context we also point out that the concept of a quasi-three dimensional speckle, 
which appears due to an extension of the autocorrelation function in the longitudinal direction of a 
transverse 2d speckle, is not applicable for the true 3d speckle, since it requires an additional space 
dimension. Then we propose a possible experimental realization for an isotropic 3d laser speckle 
potential and derive its corresponding autocorrelation function. Using a Fourier transform of that 
function, we calculate both condensate depletion and sound velocity of a Bose-Einstein condensate 
as disorder ensemble averages of such a weak laser speckle potential within a perturbative solution 
of the Gross-Pitaevskii equation. By doing so, we reproduce the expression of the normalfluid den- 
sity obtained earlier within the treatment of Landau. This physically transparent derivation shows 
that condensate particles, which are scattered by disorder, form a gas of quasiparticles which is 
responsible for the normalfluid component. 

PACS numbers: 67.85.Hj, 46.65.+g 



I. INTRODUCTION 



The study of interacting bosonic atoms in a disordered 
potential landscape, called in the literature as "dirty bo- 
son problem" [l|, has originally been introduced in the 
context of the motion of superfluid helium in porous Vy- 
cor glass [2|. Due to the frozen environment, disorder 
ensembles averages of physical observables have to be de- 
termined, vifhich depend on many system parameters as, 
for instance, the strength of a repulsive interaction be- 
tween two particles of the Bose gas as well as the strength 
and the correlation length which characterize the disorder 
potential. The main and intriguing part of the problem 
is the competition between the repulsive two-particle in- 
teraction and the localization property of disorder. From 
a theoretical point of view, the disorder potential was in- 
troduced by investigating the Anderson localization phe- 
nomenon for fermions Q. Much attention has recently 
been paid for the Anderson localization and the propaga- 
tion of bosonic matter waves in random external poten- 
tials [J]. Experimentally, the bosonic matter waves have 
been studied in the random potential produced either by 
laser speckles |5| or by an incommensurable optical lat- 
tice 0. Whereas the laser speckle disorder potential is 
created by a laser beam scattered from a diffusive glass 



plate Q , the incommensurable optical lattice is produced 
through two interfering laser beams with incommensu- 
rable wavelengths. However, one needs to remark that 
such lattices exhibit certain pathological features, which 
do not occur in genuinely random lattices, such as a tran- 
sition between localized and delocalized states, even in 
one spatial dimension [8| . In that sense the quasi-periodic 
lattices should be considered as to be quasi-random ones. 
Recent progress in different experimental realizations of 
laser speckle disorder is reported in Rcfs. [9l[l(]l|. 

According to the laser speckle theory described in the 
seminal work of Goodman [7|, lll| . the monochromatic 
light reflected from a rough surface on the scale of an op- 
tical wavelength yields many independent dephased but 
coherent wavelets which interfere at a distance, which is 
essentially larger than the wavelength. This results in 
a granular pattern of intensity that is called Gaussian 
speckle as the real and imaginary parts of the field am- 
plitude form a circular complex Gaussian distribution at 
any fixed spatial point. Details of the speckle formation 
will be considered in the next section of the paper. Here, 
we note that this distribution consists of the first-order 
statistics of the speckle disorder, while the second-order 
statistics of disorder is represented by its autocorrelation 
function. 

In order to understand the underlying physics of laser 



speckles, let us briefly describe their formation in 2d. Ob- 
ject waves are fields, which are a result of the incident po- 
larized monochromatic field reflection from a rough sur- 
face, and they are described in a plane a, /3 immediately 
adjacent to the surface in terms of a complex function 
a{a,l3) [12|- The Huygens-Fresnel principle establishes in 
the Fresncl approximation a relation between these ob- 
ject waves a{a, /3) and the complex waves A(a;, y) in the 
observation plane x, y through an integral which resem- 
bles a Fourier transformation. Hence, the wave A{x, y) is 
a result of the interference of all object waves in the x, y 
plane. As in the Fresnel approximation one assumes the 
condition z ^ (p? -I- fP'^-aiB.^I\ where z denotes the dis- 
tance between the object wave a, /? plane as well as the 
observation wave x, y plane and A denotes the light wave- 
length, the waves A{x, y) arc called to be in far field p^ . 

In the Fourier mapping of object waves for the forma- 
tion of far fields both the form and the finite size of the 
diffraction aperture A in the a, /3 plane plays a central 
role. It determines the form of the autocorrelation func- 
tion as well as its correlation length, which characterizes 
the average size of the speckle, i.e. a grain of the above 
mentioned intensity pattern. Typically, the expression 
for the autocorrelation function consists of a constant and 
a spatially varying part. The latter, which is of interest 
for various speckle applications, has one central maxi- 
mum and a set of side maxima of decaying height, which 
are separated from each other by zeros. This analyti- 
cal structure is principal in the theory of laser speckles, 
since it is the result of the Fourier transformation of the 
finite-size diffraction aperture A. Due to the existence 
of zeros, it can qualitatively not be approximated by a 
function of a Gaussian form as was assumed and even nu- 
merically derived in Refs. |13l - [la |. It is interesting that 
the experiment demonstrates an ambiguity in the follow- 
ing respect: whereas the function with zeros is exploited 
in the papers 0, [la, [l3l , the spatial autocorrelation func- 
tion is fitted by a Gaussian in Refs. p. llOl. Il8l - [2l| . Calcu- 
lating a standard deviation of the second-order moment 
of the random intensity, it was shown in Ref. [231 that 
for \d the autocorrelation function derived in Ref. j7l.[l]| 
can be well approximated by a Gaussian form. However, 
a Fourier transform of this autocorrelation function, the 
power spectral density, which is essential for the theory 
of a Bose-Einstein condensate (BEG) in an external dis- 
order potential, behaves, unlike the Gaussian function, 
as the triangle function tri(a:) = 1 — |a;| for |a;| < 1 and 
otherwise zero for any dimensionality. For Id and 2d 
this was shown by Goodman in Refs. 0, [il|, the corre- 
sponding 3d case is dealt with below in the text. This 
triangle function makes the upper limit of the integration 
in momentum space finite. For those reasons the recently 
proposed Gaussian autocorrelation function for the laser 
speckle is not suitable for a comprehensive description of 
a BEG in laser speckle disorder. 

The present paper is organized as follows. We start 
with describing the basic principles of the laser speckle 
theory in Sec. [Hi Following a scheme described in Refs. [3, 



[Tl| . we will then derive in Sec. IIIII the expressions for 
the autocorrelation function of laser speckles and their 
Fourier transforms ranging from \d to Sd with special 
emphasize on discussing both isotropic and anisotropic 
cases. The scheme of the possible experimental realiza- 
tion of the 3d isotropic speckle will be outlined in Sec. lIVI 
Note, however, that we consider in our paper a true 3d 
speckle pattern, not a quasi-thrce dimensional one of a 
transverse 2d speckle with a longitudinal depth in the au- 
tocorrelation function as described in Ref. [23l and sec- 
tion 4.4.3 of the Goodman book [il|, which has bee n ap - 
plied in many experiments (see, for instance, Ref. [22l|). 
This depth autocorrelation function concept assumes the 
existence of an additional spatial direction for the rele- 
vant speckle and can only be valid for Id or 2d speckles. 
As is further discussed in Refs. [Ill, [2J], the depth size 
is essentially larger than ones in other dimensions. Here 
we consider a 3d volume speckle with compatible speckle 
grain sizes in all spatial directions, which was already 
simulated in Refs. [13|, [ij]. Since the existing speckle 
patterns are experimentally produced mainly in a 2d ge- 
ometry, we will propose a special scheme for its possible 
realization in a 3d volume. In the subsequent Sec. [Vl the 
effect of a weak 3d isotropic speckle on various thermody- 
namic properties of a dilute Bose gas will be considered 
at zero temperature. To this end we calculate both con- 
densate depletion and sound velocity of a BEG within 
a perturbative solution of the Gross-Pitaevskii equation. 
Afterwards, in Scc. lVIi we reproduce the expression of the 
normalfluid density of a BEG in an external disorder po- 
tential obtained earlier within the treatment of Landau. 
From this rederivation we realize that condensate parti- 
cles, which are scattered by a disorder potential, form 
a gas of quasiparticles, which is responsible for the nor- 
malfluid component. Finally, we summarize and analyze 
the results obtained in the paper in Sec. I VIII 



II. FUNDAMENTALS OF LASER SPECKLE 
THEORY 

According to Refs. [3, [ll|, [I^ the circular Gaussian 
probability density function 
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for the real Ab. and imaginary Ai parts of a far-field 
A{x, y) at each point x, y with the variance 77 = a/(|^P) 
represents the background of the theory of laser speckles. 
Another basis of the theory is the A/-fold joint Gaussian 
probability density function 
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for far-fields A{x,y) at different points x^y. Here 
[Ca] is a Hermitian symmetric matrix with determi- 
nant \Ga\, whose elements are given by {CA)i,j = 



{A*{xi,yi)A{xj,yj)) for a set of far-fields [A] = 
{A{xi,yi),A{x2,y2),---A{xM,yM)} at M points of the 
x,y plane. Note that the notation (• • •) in the expres- 
sions for 77^ and {CA)i,j and throughout below in the 
text means the disorder ensemble average. Furthermore, 
one assumes that the indices i,j at {CA)i,j are taken for 
adjacent spatial positions. 

Expressions ([T|) and ^ are the result of the central 
limit theorem of probability theory [2J), which claims 
the following: if complex random variables are the sum 
of other independent complex random variables then, at 
the increase of the number of second ones, the first ones 
are distributed according to the Gaussian law. Applying 
the theorem for our case we have far fields A{x, y) as 
result of the interference of independent object waves at 
all positions of x, y plane. As we will sec below, there 
are mainly two physical conditions for the validity of the 
central limit theorem. They are related to the physics 
of providing independence of the object waves and to 
the method of their summation within the interference 
process. We will now describe both of them in more 
detail. 

A requirement for the object wave a(a, /3) to be inde- 
pendent leads to some limitations for its statistical prop- 
erties [l2l- First of all, formed after the reflection of 
monochromatic light from the rough surface, the individ- 
ual wavelet a{a, /3) should be completely polarized. Sec- 
ond, the first-order probability density of its phase should 
be uniform in the interval —n to n. And at last, the ob- 
ject wave a{a,l3) should be quasi-homogeneous, which 
means that its autocorrelation function Ca consists of a 
slowly-varying intensity la envelope and a short-range 
normalized correlation function C^: 

Ca(ai,/3i;a2,/32) = (a*(ai, /3i)a(a2, /32)) 

" V 2 ' 2 j W(a2 -ai,/32-/3i). (3) 

If we increase in this expression the range of variation 
of the correlation part, i.e. the correlation length of the 
object wave, the changing range of the intensity becomes 
smaller. However, in order to entirely satisfy the inde- 
pendence condition of the object waves, their correlation 
length in Eq. (jSj should be as short as possible, which 
means that C'^ia^P) has to be delta correlated. The lat- 
ter introduces some demands upon the properties of the 
random light scatterer, which is called in the literature 
as the diffusor. Typically, a diffusor is an optically ho- 
mogeneous transparent glass plate with no reflection cen- 
ters for light in the volume and a geometrically inhomo- 
geneous distribution of reflection centers with random 
heights on its surface. As mentioned in the introduc- 
tory section of the paper, the scattering rough surface 
generates object waves in a plane a, /3, which is closely 
situated at the surface, when the monochromatic polar- 
ized incident light transmits through the plate. Another 
realization of object waves is considered in Refs. 0, |Tl[, 
where the lateral monochromatic light was directly inci- 
dent on the rough surface. For our purpose to calculate 



the normalized speckle autocorrelation function, the op- 
tical property of the medium, from which light falls on 
the rough surface, is merely dropped from the consider- 
ation. Assuming a Gaussian probability density of the 
surface height h(a, (3) with the autocorrelation function 
Ch{oi, f3) and a variance ijf^ and assuming also a Gaussian 
probability density of object wave phases with a variance 
?7?, Goodman derived the relation j7l[ll| 



C;(a/3) = exp(-772[l-C;,(a,/3)]) 



(4) 



where 77^ = 2-Krjh/X- This function can be approxi- 
mated by a delta function (5(a_/3) when rj^ > 1 and thus 
r]h > A/2 and the mean distance between two inhomo- 
geneous h{a, /3) is larger than A. The delta functional 
autocorrelation of object waves provides their indepen- 
dence from each other. On the other hand, the Gaussian 
probability density of phases reduces to a uniform one for 
770 > 1 supporting the second requirement for the object 
waves outlined above. Therefore, the requirements for 
object waves described in the previous paragraph can be 
experimentally realized if the size of the surface inhomo- 
geneities and the distance between them are larger than 
the wave length of the light. As a next implication of the 
presented analysis we can suppose that for the outlined 
system parameters the object wave probability density 
itself may have a circular Gaussian form for the real and 
imaginary components of the wavelet a{a, /3). 

The next important step of the theory of Gaussian laser 
speckles is the formation of far fields for a given set of ob- 
ject waves. It is based on the Huygens-Frcsnel principle 
of optics, which preserves the individual wavelet picture, 
i.e., works in the limit of optics, where deviations from 
geometrical optics are small. The interference of the ob- 
ject waves yields an amplitude A{x,y), which reads in 
the above mentioned far field Fresnel approximation as 
follows 



A{x,y) = / a(a,/3)exp 
Ja 
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where we have omitted unimportant multipliers in front 
and inside of the integral. This expression resembles, 
indeed, a Fourier transformation and, thus, conserves the 
principle that each object wave contributes individually 
to the interference. 

Using Eqs. ([3]) and ([5]) it is straight-forward to derive 
the expression for the autocorrelation function 

CA{xi,yi;x2,y2) = {A*{xi,yi)A{x2,y2)) ■ (6) 
It turns out to be given by 

(7) 
for Ax = a;2 — Xi and Ay = 772 ^ J/i with 

da d/3 
(8) 
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for the novel variables a' 
and a" = a2 — cti, P" = 
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(a^a' + y/3') 
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= (ai 
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fa2)/2,/3' = (/3i+/32)/2 
/3i. Eq. ([7]) shows that 



the far fields A{x, y) are quasi-homogeneous like the ob- 
ject waves a{a,l3). This is a direct result of the Fourier 
transformation ([S]) . Another important result of this lin- 
ear transformation is the implicit proof of the above made 
supposition that object waves are circularly Gaussian dis- 
tributed. Indeed, only Gaussian distributed object waves 
can contribute through the linear mapping to Gaussian 
far fields. As we will see below the role of the so far 
uninvestigated parameter, the aperture A, will lead to 
the formation of a correlation length of the correlation 
function C^(a;, y). 

The autocorrelation function Ci{x,y) of the far-field 
intensity I{x, y) — \A{x, y)p can be calculated using the 
Wick theorem for variables distributed according to the 
Gaussian law. A simple calculation gives the expression 



Ci{x,y)^{I)^l + \C'Ax,yW 



(10) 



where the normalized autocorrelation function for the 
far-field C\{x,y) is defined as C^(x-, j/)/C^(0, 0), where 

As already mentioned above, if in Eq. ([3]) the autocor- 
relation function of object waves C^(a^/3) is delta corre- 
lated, then the intensity function of these waves /a (a, (3) 
can be approximated as a constant. Assuming that the 
a,/3 plane is close to the rough surface, one can write 
|a(a,/3)| = K|P(a,/3)|, where P(a,/3) is the incident to 
the glass plate light wave and k is the average refiectiv- 
ity of surface, for each position a, 13. Then the intensity 
of the object waves at a, (3 is determined by the relation 
Ia{a, f3) = K^|P(a,/3)p. Therefore, the expression for 
the normalized far-field autocorrelation function reads 



c'Aix,y) 
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III. SPECKLE AUTOCORRELATION 

FUNCTION FOR APERTURES IN Id TO 3d 

DIMENSIONS 

As already mentioned in the introductory section, the 
investigation of a BEG in the laser speckle disorder has 
found much attention from both a theoretical and an 
experimental point of view. In particular, a variety of 
isotropic and anisotropic speckles have been the subject 
of these works. Motivated by this interest, we will de- 
scribe in the present section the derivation of the speckle 
autocorrelation function for different apertures ranging 



from one to three dimensions by generalizing the appro- 
priate expressions from the previous section to these di- 



mensions. 



A. Real space 

Due to the analytic form of Eq. pTj) . we can take the 
intensity of the incident wave |P(a, /3)p to be unity over 
the whole aperture region of the a, /3 plane. Writing the 
function |Pp in the form |P|^^, where d is the space 
of dimensionality and A is the form of the aperture, we 
have the following expressions: 



l^(«,/3)|L. 



rect I -— 



^ 

Lp 



rcct {^\ (12) 



for the 2d anisotropic rectangular aperture with sizes La 
and Lp, where the function rect(a;) = 1 for |a;| < 1/2 and 
zero otherwise; retaining in Eqs. (jl2[) only the first rect(a;) 
function and equating La = L one obtains the expression 
of |P(a)lidinv for the Id interval aperture of the size L; 
the analytic form of \P{a,(3)\2d,q_At for the 2d quadratic 
aperture of the size L is obtained from Eqs. ((T2|) if we 
specialize this equation according to La ^ Lp = L; 



l^(«,/3)lL, 



(13) 



for the 2d isotropic circular aperture with the diameter 
D and r = \J a^ -\- /3^, where the function circ(a;) = 1 for 
Ixl < 1 and zero otherwise; 



P_ 



rect 



l-P(":/5:7)l3dacpi = rect ( — ) rect 

'^(14) 

for the 3(i anisotropic rectangular parallelepiped aperture 
of sizes La, Lp and i^; the expression |J^(a,/?,7)l3rf cub 
of the 3d cubic aperture of the size L is obtained from 
Eq. (dH) by setting La = Lp 



Ij^y - — ■ 1j'^ 



I^(«,/3,7)IL. 



sph 



(15) 



(11) for the 3d isotropic sphere aperture with the diameter D 



and r = \l o? + /3^ + 7 



2. 
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(16) 



for the 3d anisotropic cylinder aperture with the diam- 
eter of circle D, r = \/ oP- + /3^ and size L^ along the 7 
axis. Substituting the expressions ([T2|) -(fT6 |) of the |P|^ ^ 
function in Eq. (|11[) and calculating the respective inte- 
grals, we obtain the corresponding expressions for the 



correlation function |C^p: 



\C'A^x,^y)\li,, 
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Xz 



(17) 



where sinc(?/) = sin(7ry)/(7ry), for the 2d anisotropic rect- 
angular aperture with z being the distance between ob- 
ject wave and far field planes; retaining in this equa- 



tion only first 
Aa; and assuming 



sine (y) function, the dependence on 
La = L one obtains the expression 



C^(Aa;)|^j^ jjj^ for the Id interval aperture with z being 



the distance between object wave and far field intervals; 
the expression \C'ji^{/S.x, /S.y)W^ ^^ for the 2d quadratic 
aperture one can derive from Eq. (J17p for La ~ Lp ~ L; 
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f TiDr 



TiDr 
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(18) 



where Ji (x) is a Bcssel function of the first kind and of 
the first order, for the 2d isotropic circular aperture with 
r = V(Aa;)2 + (Ay)2; 



|C;i(Ax,Ay,Az)|2^,,,pi 



La^X 

Xz 
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for the 3d anisotropic rectangular parallelepiped aperture 
with z as the distance between the object wave and the 
far field volumes; the expression \C'^{Ax,Ay,Az)\'^j^^^^ 
for the 3d cubic aperture is obtained from Eq. ([T^ by 
assuming La ~ Lp = L^ ~ L; 
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for the 3d isotropic sphere aperture with 
V(Ax)2 + (A2/)2 + (Az)2; 
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(21) 



for the 3d anisotropic cylinder aperture with 
V(Ax)2 + [Ayy. 

Expressions of the autocorrelation function for a 2d 
quadratic aperture |C^(Aa;, Aj/)|2^ ^^ and for a 2d cir- 
cular aperture in Eq. (jlSp are derived by Goodman in 
Refs. P, [ll[. As can be seen from the formulas for other 
cases of the aperture, they are closely related to both 
of these Goodman cases of the aperture. However, the 
derivation of the 3d isotropic sphere autocorrelation func- 
tion ()20p . which is a result of the present paper, required 
some additional effort. 

The analytical forms of the autocorrelation functions 
|C^P are similar in every spatial direction. They have 
one central maximum and a set of side maxima of decay- 
ing height, which arc separated from each other by zeros. 



As was pointed out in the introductory section, it is ob- 
vious that these forms can not be fitted by a Gaussian. 
The argument of the autocorrelation function, which cor- 
responds to its first zero, provides the correlation length 
of the disorder, i.e. the average size of the speckle grain, 
for the appropriate spatial direction. Denoting it by 5x 
we have, for instance, for Id speckle 



(. Xz 
ox — — . 
Ij 



(22) 



The main interest of the present paper is the 3d spher- 
ical aperture of Eq. ([20| since we will carry out the 
calculation of BEG properties for this particular case 
of laser speckles. Numerically solving the equation 
sin(a;) — x cos(a;) = we find first its solution to be at 
Xc — 4.493; thus the disorder correlation length is given 
by re = 1.4302 Xz/D. 

In order to establish a physical meaning of 5x we intro- 
duce the "wave number" fcoff, which is related to the vec- 
tor a, j3 in the above Fourier transform formulas, by the 
relation k^s = 271x1 (Xz). If we substitute in it 5x from 
Eq. p2|) then we obtain fccff = 2-k jL. For a circular and 
a spherical aperture the "wave number" is fcoff = 27r/Z?. 
However, the sense of fcoff is in an uncertainty of the wave 
vector when the problem of wave propagation is solved 
in the restricted area. It is well known that in this area 
the wave vector is determined within the resolution fcoff . 
Therefore, we can say that the origin of a speckle grain 
with a correlation length bx as its size represents the spa- 
tial uncertainty in the determination of far fields, which 
is introduced by the finite size of the aperture. 



B. Fourier space 

For many applications the Fourier transform of the 
far-field intensity autocorrelation function, or the power 
spectral density, of the speckle is of considerable inter- 
est. In the literature on laser speckle theory [3, [ill it is 
defined according to 



C/(k) = / C/(x)e 



-i27rkx 



d^x. 



(23) 



Substituting in it Eq. pU)) for C/(x) one obtains 

C,(k) = (/)2[<5(k) + |C;i(k)p] . (24) 

In the perturbative considerations of BEG in the speckle 
potential the Fourier transform | C^ (k) ^ plays the central 
role. It has the following expressions for the real space 
autocorrelation functions taken from Eqs. pT)) - (l211) : 
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where the triangle function is defined as tri(.T) = 1 — \x 
for \x\ < 1 and zero otherwise, for the 2d anisotropic rect 
angular aperture; the expression |C^(k)| 



L,i„v for the Id 



interval aperture one can get from Eq. (1251) by assuming 
Lq, = Li3=L,kx = ky = k and taking the square root of 
its right-hand side; the expression |C^(k)|2(j^qdt for the 
2d quadratic aperture is obtained from Eq. ([25]) with the 
assumption L^ = Lp = L; 



|C'A(k)|2d,crc 




(26) 



for the 2d isotropic circular aperture with k = i//s^. + k^: 
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for the 3d anisotropic rectangular parallelepiped aper- 
ture; the expression |C'^(k)|3^ ^^i, for the 3d cubic aper- 
ture is obtained from Eq. (|27p by specializing La ^ Lp = 

Ly = L'^ 



_ 3 f2\z 



|C;i(k)|^,,p, = 3(^) (6^-126+16) (28) 



for the 3(i isotropic sphere aperture with h — 2k\z/ D 
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for the 3d anisotropic cylinder aperture with 
Jkl + kl 



k = 



The equation for the quadratic aperture IC'aC^)!^ odt 
and Eq. ([^5)) have been derived by Goodman in Refs. [7], 
\Vi\ . Other expressions of |C^(k)p, except for the 3d 
isotropic sphere aperture case, can be obtained by us- 
ing these formulas. Eq. ([28|) is a result of this paper. 
In all our formulas for the anisotropic aperture we have 
assumed that the size deviation of the aperture with re- 
spect to its average isotropic size is essentially less than 
the distance z. 

As is seen from the formulas of |C^(k)p expressed 
through the triangle function their value becomes zero 
when their argument is unity. For the 2d circle and the 3d 
sphere apertures |C^(k)p is zero for kXz/D — 1. Hence, 
the wave vector of the Fourier transform autocorrelation 
function only varies in a finite interval from zero, in con- 
trast to the case for a Gaussian function. This fact is 
another reason why the speckle autocorrelation function 
can not be approximated by a Gaussian form. 



It is worth to discuss the expression |C^(7')| = 
sinclkLr) with fc f, ~ D/Xz for the autocorrelation func- 
tion used in Ref. [231 fo^ ^^^ 3d isotropic aperture. It is 
similar to our correlation function |C^(Aa;)|?^ jj^^ for the 
Id interval aperture. The authors of Ref. [25[ claim that 
this expression is vahd for z ~ (a^ + /3^)max/A, which is 
outside of the far- field limit. However, that limit destroys 
the fundamentals of the Gaussian speckle theory as they 
are described in Sec. [Ill Therefore, it is unclear whether 
|C;i(r)|2 of Ref. 



25[ is related to laser speckles or not. 



Furthermore, we discuss the definition of the speckle 
correlation length to be the width at the half value of the 
maximum of |C^(r)p for r = 0, when the last one is ap- 
proximated by a Gaussian function. Probably, this defi- 
nition was introduced first by Modugno in Ref. [26| , when 
he considered |C'A(^^)lid inv ^^ ^^^ found in Ref. [2y| 
that the correlation length is given by Sx = 0.88Xz/L, 
while from Eq. (P^ the exact value turns out to be 
Sx = Xz/L. It is interesting that the Gaussian |C^(r)p 
has been obtained in the numerical simulation of the 3d 
isotropic laser speckle in Refs. [Ij, [ij which should be 



compared with the exact |C^ 



WiL.sph inEq. 



(HOD, with 



the correlation length re = l.lXz/D, however, the exact 
one is Tc = 1.4302A2/Z?, see the discussion after Eq. (|22p . 
It seems that we can explain the reason why the authors 
of Refs. [Ij, 13 obtained the Gaussian form of |C^(r)p. 
They used the speckle simulation method proposed by 
Huntley in Ref. [27| which we briefly review for the 2d 
case. Let us consider to this end two square planes a, /3 
and X, y with the same size L. According to the Hunt- 
ley method one uses Eq. (O in order to perform a double 
Fourier transformation. In the first inverse Fourier trans- 
formation the complex object waves a{a, j3) on the mesh 
points in the a, /3 plane are simulated through the given 
Gaussian distributed complex random waves A{x, y) on 
the mesh points in the x, y plane. Afterwards, one cuts 
by a circle with radius D/2 the a, (3 region of the ob- 
tained a(a, 13) such that it vanishes outside of this region. 
In the second direct Fourier transformation the derived 
complex waves a(a, 13) form the final complex far-fields 
A{x,y). Huntley has investigated in Ref. [23| only the 
first-order statistical property of the simulated pattern, 
i.e. the probability density of the intensity, and showed 
that it corresponds to the theoretical laser speckles of 
Ref. [3|. However, the proposed simulation method can 
drastically deviate in the second-order statistical prop- 
erty of a speckle, i.e. its autocorrelation function, from 
the theoretical one. 

Indeed, in accordance with the theory of a speckle au- 
tocorrelation function as presented in this section, after 
the first Fourier mapping the object waves a(a, j3) ac- 
quire a correlation with the correlation length 5a = Sp = 
Xz/L, where L is size of the square x, y plane. More 
precisely, now the function C^(a,/3) is not delta corre- 
lated. However, according to Eq. ([3]), the broadening of 
the C^(a,/3) function reduces to a changing of a con- 
stant character of the Ia(oij (3) function to one of varying 
in space in the a^(3 plane. Substituting this function of 



Ia{a,l3) in Eq. (J9]) and integrating over a and /? gives 
the function C'j^{x,y) which may quaUtatively be differ- 
ent from the one discussed in this section. 

A simulation method, which is consistent with the 
above laser speckle theory, is described in the book of 
Goodman lllll . There are other numerical methods in 
Refs. [2a [23, in which the exact form of the real space 
autocorrelation function is used to generate the speckle 
pattern. In particular, one of such methods was exploited 
for the simulation of Id speckle in Ref. [30l |. 



IV. EXPERIMENTAL REALIZATION OF 3d 
ISOTROPIC SPECKLE 




As was already mentioned in the introductory section, 
we consider here a true Sd speckle, not the quasi-threc di- 
mensional one consisting of a transverse 2d speckle with 
a longitudinal depth in the autocorrelation function as 
described in details in Ref. [23| and section 4.4.3 of the 
Goodman book [ll| and applied in many experiments. At 
a first glance, it seems exotic and unrealistic to experi- 
mentally realize such a Sd volume speckle pattern. How- 
ever, in the present section we will describe the physical 
principle how it can be generated. 

In the typical 2d geometry of the experimental realiza- 
tion of a speckle a lens, which collects the incident light, 
is installed close to the glass plate such that its focal 
plane coincides with the far- field plane |22|. This idea 
of a speckle formation in the focal plane can be gener- 
alized to a full 3d geometry, when the speckle is formed 
in the focal point, i.e. the focus, of an empty ellipsoidal 
optic cavity according to the scheme displayed in Fig. 1. 
Let us consider that cavity, whose inside surface refiects 
absolutely the light emitted from one of its focus (point 
A) and collects it at the second focus (point B). Two 
laser beams (thick yellow arrows) are incident through 
holes in the cavity surface into the small metallic sphere, 
i.e. the reflector, with absolute light reflection, located 
in the center of the glass sphere A. It is assumed that 
laser beams cover the entire surface of this reflector and 
a BEG is deposited in the sphere B, which is located at 
the second focus of the ellipsoid. 

The glass sphere A additionally contains the located 
randomly light scattering centers, for instance, absolutely 
light reflective metallic polyhedrons with a random aver- 
age size of each facet. The theory how to derive the 2d 
object wave autocorrelation function C^(a,/3), described 
in Refs. [3,|li|, can be easily generalized to the derivation 
of (7^(0,/?, 7) for such a 3d case with the same expres- 
sion Q. However, now this expression is a function of 
3d other quantities. The condition, at which C^(a,/3,7) 
becomes delta correlated and the object waves are inde- 
pendent, is the same as for C^(a,/3). Therefore, if the 
mean distance between these light scattering centers and 
the average size of polyhedrons are larger than the light 
wavelength, then C^(a, /3, 7) will be delta correlated. On 
the other hand, each sphere with radius (a^ -f /3^ -f-7^)^/^ 



FIG. 1: Cross section of the ellipsoidal reflective cavity with 
spheres A and B in its focuses. Focus A contains a small 
size absolute spherical light reflector in the center and vol- 
ume optic inhomogeneities. Incident laser beams (thick yellow 
arrows), after reflection from the reflector, scatter addition- 
ally from inhomogeneities producing individual wavelets (thin 
yellow lines), which are collected in focus B, where a BEC is 
deposited. 



and with the same center as the sphere A will be the ob- 
ject wave volume, whereas for comparison for 2d we had 
'A a, P object wave plane. 

Incident laser beams, after reflection from the reflector, 
scatter additionally from scattering centers and produce 
individual and independent wavelets, the object waves 
a{a^ 13, 7), indicated via thin yellow lines in Fig. [l] which 
are collected in the sphere B, where a BEC is deposited. 
For the presented geometry the far-field condition is sat- 
isfied, since a distance z between the object wave and the 
far-field volumes, i.e. the length of each wavelet trajec- 
tory between two focuses of the ellipsoid, is larger than 
the size of the object wave sphere A. 

The described scheme can be generalized for the ex- 
perimental realization of any 3d anisotropic speckle. To 
this end one only needs to change the form of the spher- 
ical aperture A, which contains the glass and the light 
scattering centers, into a suitable one listed in the previ- 
ous section. The spherical form of the metallic refiector 
retains unchanged. 

At the end of this section, it is worthwhile to discuss 
the possible realization of a 3d volume speckle pattern 
using 2d plane speckles. Such a scenario presumes a 3d 
speckle as a result of the sum or, more clearly, as a linear 
interference of two and more 2d speckles. While the- 
oretically this scenario is discussed by Pilati at al. in 
Ref. [l4| . the experiment, in which two perpendicular 2d 
speckle planes form a 3d speckle pattern, was realized 
in Ref. [3l| by Jendrzejcwski et al. Instantly the ques- 
tion arises whether the random pattern realized in such 
a way belongs to the class of speckles or not. In spite 



of an additional theoretical analysis, which is required to 
answer that question in detail, the following argument 
shows that the possible conclusion is negative. Indeed, 
according to the fundamentals of the laser speckle the- 
ory of Goodman, Refs. 0, [ll|, and Dainty, Ref. 12|, see 
also Sees, im and mil of this paper, the correlated speckle 
pattern in any dimension, except the one described in 
Ref. [23| and its analogue for Id (see next paragraph). 



is a result of the Fourier transform over the restricted 
aperture object wave region in the same dimension. This 
means that a 3d volume speckle can be obtained only 
by a 3d object wave volume. Physically it means that 
the single connected spatial domain of each 3d speckle 
grain, which is a result of 3d correlations, can not be 
obtained by a linear combination of randomly sized and 
independent 2d speckle grains. By that reason, the true 
3d speckle cannot be obtained even by a combination of 
quasi-three dimensional speckles, which we discussed at 
the beginning of this section. 



V. BEC DEPLETION AND SOUND VELOCITY 
IN WEAK 3d ISOTROPIC SPECKLE 

The interaction potential of light with an atom at po- 
sition r is determined by the far-field intensity /(r) ~ 
\A(r)'^\ and has the form V{r) = tl{r), see for instance 
Refs. [22, [23, where the constant i is a function of the 
atomic and light characteristics. At the derivation of 
V{r) it was assumed that the incident laser wave does 
not induce an atomic electron interlevel transition, but 
merely deforms the atomic ground state. 

It is convenient to define the interaction potential as 
V{r) = Vo+ AV"(r), where Ay(r) = V(r) - Vq and 
Vb = (/). Using the obvious property (AT^(r)) = 0, a 
simple calculation shows that 



{V{r')V{r'+v))=V,' 



1 



(Ay(r')AF(r' + r)) 



^0 



(30) 



and, therefore, we have the following relationships be- 
tween the laser speckle autocorrelation and the disorder 
potential correlation functions: 



|CKr)P = {V{r')V{r' + r)), 

2 _ (AF(r')Ay(r' + r)) 



\C'A{r)\' 



^0 



(31) 



Our interest is a Bose gas with a contact interaction. 
Taking into account that, according to the novel defini- 
tion of V{r), the chemical potential for the ground state 
of BEC will be renormalized according to /x — > /i — Vq, 
the Gross-Pitaevskii equation (GPE) reads 



2m 



V2 + Ay(r)+5|*(r)|2 



*(r) = 0. (32) 



Here g — inh^a/m denotes the strength of the contact 
interaction with the scattering length a. 



Under the assumption that the disorder potential is 
weak, one can expand the solution 



*(r) = Vo + V'i(r)+V'2(r) 



(33) 



and solve the GPE ([5^ perturbatively in the respective 
order of Ay(r) [34|. For the ground state all functions 
of the expansion as well as ^(r) are real. In this way 
the problem is reduced to find the total particle density 
n ~ (\I'(r)^) and the condensate density no = (^(r))^. In 
particular, the lowest order expression for the condensate 
depletion reads 



72 — riQ 



no 



d^k 



i?(k) 



(27r)3 [n^k^/2m + 2ng] 



(34) 



where we introduced the literature notation i?(k) = 
In order to further apply our formula Eq. ((28)) for the 



3d isotropic autocorrelation function |C^(k)| 



3rf,sph' 



one 



needs to make a remark. According to the definition in 
Eq. (P5)) , the Fourier transforms of autocorrelation func- 
tions carry a physical dimension. In particular, the cor- 
relation function |C'^(k)|3^ ,, j^, calculated with Eq. (^5]) . 
is proportional to the inverse volume of the 3d isotropic 
aperture 3/(47r)(2/I?)'^ times (Az)'^. If we introduce the 
correlation length as ct = \z/ D, then the proportionality 
factor is 3{2aY /{Att). In the following, we assume that 
|C^(k)|3^ g jj is already normalized by that factor. 

It is convenient to introduce also the BEC co- 
herence length according to ^ = [It? /{2mng)Y/'^ ~ 
l/voTma. Substituting the normalized correlation 
function i?|C;i(k)|2^^pi^ from Eq. ^ in Eq. ^ 
and performing the integration, we get the expression 
n — no = nnufiu/S^), where the depletion n-^u = 
[m'^R/{8n^^^h )] \frifa was obtained by Huang and Meng 
in Ref. [33| (see also Ref. |3J]) for delta correlated disor- 
der i?(r) and the condensate depletion function is defined 
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^ arctan . ^ 
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(35) 



The function f{(j/^)^ which is depicted in Fig. [21 has the 
following asymptotics for small a /S, 



f 



1- 



14x/2 
37r 



I8V2 /o- 



and, correspondingly, for large a/^ 

3 



/(f 



1 



25/2^ 



- r- 

3 \a 



10 Vcr 



+ 



(36) 



(37) 



Introducing the appropriate correlation length for each 
aperture, as described in Sec. IIIIl one can show that. 





1.0 



FIG. 2: Condensate depletion function /(o-/^) from Eq. (|35 



FIG. 3: Sound velocity function s((t/^) from Eq. (f39l 



when this correlation length tends to zero, then the cor- 
responding correlation function |C^(r)p tends to the 
delta function. 



for small ct/^ 



|CA(i')lL,sph in the limit a ^ 0. 



The same behavior has our function 
Therefore, vife should 
reproduce the Huang and Meng result tt-hm for the con- 
densate depletion in this limit. Indeed, when cr/^ — >■ 
we read off from Eq. (|55|) that one obtains /(cr/^) -^ 1. 

For the Sd isotropic Bose gas with contact interaction 
the normalfluid density njq is determined by the equation 
TIN = 4(n — ri,o)/3 (see Sec. IVII below and Ref. [32| as well 
as the references therein) , from which ntq is proportional 
to the function of f{<7/S,). 

In Ref. [32| the sound velocity of a dipolar BEC in a 



weak external disorder potential is calculated within a 
hydrodynamic approach. To this end a general deriva- 
tion was performed which is applicable for an arbitrary 
interaction potential. For an isotropic id system with 
contact interaction it has the form; 



c 
Co 



c/3fc 



i?(k) 



(27r)3 (a^k2/27 
(?i^k2/2m + 2ng) 



+ 2ngf- 



(38) 



where Cq = {ng/mY''^ is the sound velocity in a sys- 
tem without disorder and the scalar product between the 
sound direction q and the direction of wave propagation 
k has the form qk = cos •& for an isotropic system. 

Calculating the integral in Eq. (|55)) . we obtain c/c{) ~ 
1 + riHMs(o'/C)/(2ft), where the sound velocity function 
reads 



23/2 g 



14 

y 



28^2 \, / e 



' ^.arctan^ ^ 



(39) 



3V2cr "''^""" VV2cr, 
It is depicted in Fig. [3] and has the following asymptotics 



<j\ 5 23/23 fa 



and for large cr/^ 






3 \a 



23/262 /a 
9^ U 



18 Icr 



(40) 



(41) 



respectively. 

Again, when the correlation length cr — > and thus the 
correlation function IC'aC'')!^ sph ^^ delta correlated, we 
reproduce the result s(a/^) ~ 5/3 of Ref. [33, obtained 
for delta correlated R{r). 

As shown in this section, the finite range of integration 
for the vector k, when a Fourier transform of a speckle 
correlation function is being applied, essentially simpli- 
fies the analytic calculation of the BEC properties. This 
is an essential advantage of applying the laser speckle 
theory to the BEC investigation. Conversely, due to the 
infinite limit of integration on k, the Gaussian disorder 
correlation function, which is often used in the literature, 
introduces some difficulties in its application to the BEC 
theory. 



VI. LANDAU DERIVATION OF 
NORMALFLUID DENSITY 

Let Kq be a reference frame, and K a second frame 
with relative velocity — v with respect to K^. According 
to the Galilean transformation in classical mechanics, the 
energy Eq of a system in the frame Kq and its energy E 
in the frame K are related to each other by: 



M ., 
E = Eo-Po^ + Y^ 



(42) 



where Pq and M are the total momentum and the mass 
of the system, respectively. 
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Following to Refs. (36|, |37[ let us assume that, at tem- 
perature T = 0, the condensate is in rest, i.e., in the 
frame Kq^ and its energy is _Eo = with momentum 
Po = 0. If one quasiparticle with mass m appears in 
the condensate with energy e(p), where p is a momen- 
tum of the quasiparticle, then in the frame Kq energy 
and momentum now become Eq — e(p) and Pq = P- 
Hence, from Eq. (|1^ the energy E in frame K will be 
E ~ e(p) — pv + A/v^/2 and the energy of the quasi- 
particle in frame K after a Galilean transformation has 
a form e(p) — pv. 

According to the Landau two-fluid theory [3^, [331 of 
liquid helium II, a gas of quasiparticles, for instance 
phonons, constitutes the normalfluid density at low tem- 
peratures. For r = no quasiparticles exist, thus the 
helium is entirely supcrfluid. If a gas of quasiparticles ap- 
pears in the system for finite but low temperatures, which 
has zero center mass velocity in the frame Kq and moves 
with constant velocity — v with respect to the frame K, 
in which the helium liquid is in the rest, then the total 
momentum of the gas per volume in the frame K is given 
by 



^ = /piV(s(p)-pv)^ 



(43) 



where A^(e(p)) is the average occupation number of states 
by phonons with energy £(p). Eq. (|43l) describes the ther- 
modynamic property of a gas of phonons. However, it can 
be generalized to our BEG system in the external disor- 
der potential at T = 0, if we assume that, after scattering 
with the disorder, particles of the condensate become the 
quasiparticles of the normalfluid density. It is clear that 
it occurs when the disorder is attached to the frame Kq- 
To this end we replace in Eq. (|43|) the thermodynamic 
quantity A^(e(p) — pv) by the quantum one |^(p— tov)^, 
where the wave function is a solution of the GPE with 
disorder and written in momentum representation. After 
that we average both sides of Eq. (|43|) over the disorder 
ensemble. The obtained mean square of the modulo of 
the wave function is now homogeneous in space, so it can 
be expressed in terms of the energy of a quasiparticle, 
as the Hamiltonian is commutative with the momentum 
operator and thus the eigenfunction of the latter can be 
taken as the eigenfunction of the former [33| . Recalling 
that the expression for the total density is n = (\I'^), we 
obtain 



(P) 
V 



pn(e(p) -pv) 






(44) 



In order to derive the expression for the normalfluid 
density we expand the integrand of Eq. (|44p in power of 
pv and, in the limit v — > 0, retain only its first two terms. 
After integrating over the directions of the vector p the 
zeroth order term of this expansion disappears. Thus one 
obtains 



V = - /p(pv) 



^dn{e{p)) cP'p 
' <i£(p) (27r;i)3 



(45) 



This expression is the main result of the normalfluid den- 
sity Landau theory, when the two replacements (P) by 
P and n(e(p)) by A''(e(p)) arc performed. 

Taking into account that p(pv) = p^^, the expression 
for the normalfluid density reduces to 



Pn = 



Pz 



dn{e{p)) d^p 
rfe(p) {2Trh)3 



(46) 



From Eq. (|M)) wc have the expression of the total density 
Fourier transform 



n(£(p)) = (27r)3no<5(k) 



noRCk) 



(n^k2/2m + 2ngy 



(47) 



in first order of R(k), from which the energy of the 
quasiparticles follows to be £r(p) = p^ /2m + 2ng, where 
p = Kk. Substituting e(p) in Eq. (|46p and performing 
its integral by parts and using in the obtained expression 
n{e{p)) from Eq. (|47p . one gets 



Po 



d^k 



plR{k) 



{2T:)^p^{h\^/2m + 2ngy 



(48) 



where po = mno. 

It is interesting that there is the relationship £b(p) = 
£^/2(p)p/(2m)^/2 between our e(p) and the Bogoliubov 
quasiparticle energy ^^(p). If we use this relation, then 
we obtain 



Pn = -r 



P_a 
4 



d^k p^plRik) 



(27r)3 



9 4 



(p) 



(49) 



This expression without the prefactor 1/4 coincides with 
Eq. (19) of Ref. [35| for the normalfluid density Pn^LR, ob- 
tained within the linear response approach, if wc replace 
V J d'^ k / {2tt)^ by J^w- The prefactor 1/4 appears from 
the relation between e(p) and EbIp)- For a 3d isotropic 
BEG system we have p^ = p^ ~ pi and p^ = 3p1 . Multi- 
plying the right-hand side of Eq. (j48|) with 3 and cancel- 
ing 'Spl and p2 in the numerator and the denominator, 
we obtain Eq. (p4|) . therefore, n — no = ipuXRli'^rn) j35l |. 

It is worth to discuss the validity to use the Landau ap- 
proach for BEG with the disorder. According to a remark 
in the text book [33] the Landau approach should not be 
applicable for such a system. Indeed, the applied Lan- 
dau derivation of the normalfluid density presumes the 
validity of the quasiparticle concept (see, for instance, 
Refs. [3a [331), in which there are no collisions not only 
between quasiparticles but also of last ones with the ex- 
ternal disorder potential. More exactly, according to this 
concept quasiparticles should be well deflned and their 
gas should be ideal. 

In our case, effective quasiparticles with the mean-fleld 
energy £(p) and the quantum state distribution at tem- 
perature T = 0, represented by the total density n(e(p)), 
appear in the system after the disorder ensemble average. 
However, after this averaging the real space is homoge- 
neous and there is no reason for the gas of effective quasi- 
particles to be not ideal. Hence, if for the conventional 
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quasiparticles the source of their appearance is the low 
temperature, here it is the scattering of the condensate 
particles with the disorder and then their excitation and 
departure from the condensate. This physical conclusion 
naturally arises from the Landau derivation of the nor- 
malfluid density. 



VII. SUMMARY AND CONCLUSION 



ity of a BEC in a weak speckle disorder within a per- 
turbativc solution of the Gross-Pitaevskii equation. At 
the end, we reproduced the expression of the normalfiuid 
density obtained earlier within the treatment of Landau. 
This physically transparent derivation showed that con- 
densate particles, which are scattered by disorder, form 
a gas of quasiparticles which is responsible for the nor- 
malfiuid component. We have justified the validity of the 
Landau approach to our BEC system with disorder. 



At first, we have summarized the derivation of the au- 
tocorrelation function of the laser speckle in Id and 2d 
following the seminal work of Goodman. We showed that 
a Gaussian approximation of this function, proposed in 
some recent papers, is inconsistent with the background 
of laser speckle theory. Then we have proposed a possible 
experimental realization for an isotropic 3d laser speckle 
potential and derived its corresponding autocorrelation 
function. Using a Fourier transform of that function, we 
calculated both condensate depletion and sound veloc- 
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